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Abstract. We consider a curved chain of nonlinear oscillators and show that the 
interplay of curvature and nonlinearity leads to a number of qualitative effects. In 
particular, the energy of nonlinear localized excitations centered on the bending 
decreases when curvature increases, i.e. bending manifests itself as a trap for 
excitations. Moreover, the potential of this trap is double-well, thus leading to a 
symmetry breaking phenomenon: a symmetric stationary state may become unstable 
and transform into an energetically favorable asymmetric stationary state. The 
essentials of symmetry breaking are examined analytically for a simplified model. We 
also demonstrate a threshold character of the scattering process, i.e. transmission, 
trapping, or reflection of the moving nonlinear excitation passing through the bending. 



1. Introduction 

Recently, determination of the dynamical properties of physical systems with nontrivial 
geometry has attracted a growing interest because of their wide applicability in various 
physical and biophysical problems. Examples are biological macromolecules, nanotubes, 
microtubules, vesicles, electronic and photonic wave-guide structures [0], §, [| |], [5], || ; the 
gene transcription is usually accompanied by a local DNA bending |7|, |[ || ; the electronic 



properties of carbon nanotubes drastically depend on their chirality [10]; T-shape 
junctions were recently proposed []ll] as nanoscale metal-semiconductor-metal contact 
devices; periodically spaced, curved optical waveguides were proposed for observation 
of optical Bloch oscillations pi. 



The geometry manifests itself in particular in creation of linear quasi-bound states 



(see e.g. fl3| , |T4]]). On the other hand, it is well known that the balance between 
nonlinearity and dispersion provides an existence of spatially localized soliton-like 
excitations. Taken together these two localization mechanisms compete and one may 
expect that the interplay of geometry and nonlinearity can lead to qualitatively new 
effects. 

An example of nonlinearity-induced change of the geometry of the system can be 
found in Ref. |15| where the classical Heisenberg model on an infinite elastic cylinder was 
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considered. It was shown that a periodic topological spin soliton induces a periodic pinch 
of the cylindrical manifold. These results are applicable to microtubules and vesicles 
comprised of magnetic organic materials [16]. The stationary lasing of the modes in a 



microdisk was investigated in Ref. []T7| where it was shown that nonlinear whispering 
gallery modes exist in the microdisk above some finite pumping thresholds. 

In this paper we consider a curved chain of nonlinear oscillators which is described 
by the nonlinear Schrodinger (NLS) model with the coupling coefficients depending 
on the distance between oscillators. This model has been in particular used as an 
approximation to the nonlinear dynamics of biological macromolecules, such as proteins 
|T8| , [19, [20| and DNA [21]. It was shown that nonlinear excitations should play a 



key part in the biological functioning of such molecular chains. However, usually they 
are investigated on the assumption that the chain has straight rod-like configuration. 
On the other hand, thermodynamical properties of these macromolecules cannot be 
understood without taking into account their flexibility and existence of non-trivial 
spatial configurations. Moreover, there is increasing evidence that non-trivial spatial 
configurations should play a part in nonlinear dynamics of the macromolecules. In 
particular, the DNA molecule is usually hooked [0, |8|, || by the RNA-polymerase which 
binds to the promoter segment of DNA and then initiates the DNA transcription. 

This raises the interesting question of whether a non-trivial spatial configuration of 
the chain could introduce some new qualitative effects into the nonlinear dynamics of the 
system? In the present paper we demonstrate some of such effects on the example of the 
chain with a single parabolic bending. In Sec. 2 we describe the model and present the 
results of numerical calculations of stationary excitations. In particular, we discuss the 
multistability phenomenon and curvature-induced symmetry breaking in the spectrum 
of nonlinear excitations. In Sec. 3 we analyze analytically a minimal model which gives 
qualitatively the same behavior as the original model. In Sec. 4 we consider dynamical 
properties of nonlinear excitations, in particular scattering of the moving excitations by 
the chain bending. Section 5 presents the concluding discussion. 

2. System and equations of motion 

The model we study is described by the Hamiltonian 

H = Y,U\*Pn\ 2 - E Jn, m r n ^m-h^nA , (1) 
n y m^n J 

where ip n (t) is the complex amplitude at the site n (n = 0, ±1, ±2, . . . ). For instance 



in the case of the two-strand DNA model ||22|| , ip is the complex amplitude of the 
base-pair stretching vibration. We assume that the excitation transfer J n m in the 
dispersive term depends on the distance in the embedding space between the sites n 
and m: J n , m = J(\f n — f m \), where the radius- vector r n = (x n ,y n ,z n ) characterizes 
spatial position of the site n. Of course, when J n m models electromagnetic coupling 
between dipoles, this coefficient will depend also on the relative angles. However, in 
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Figure 1. The curved chain under consideration and the corresponding effective 
double-well potential U and mass m, created by the bending. 



some cases such dependence is not present. E.g., for the semiflexible ribbon polymers 
23] the dipoles are orthogonal to the plane in which the chain can be curved, so that 



the corresponding interaction would depend only on the distance between the dipoles. 
From the Hamiltonian (|l|) we obtain the equation of motion in the form 

i-T.^n - tyn + X] J n,m1Pm + WVn = ■ ( 2 ) 
rn=£n 

The Hamiltonian H and the number of excitations (quanta) N = J2 n iV'nl 2 are both 
conserved quantities. The sites are assumed to be equidistantly placed on a planar 
envelope curve. 

We consider the properties of nonlinear excitations in the vicinity of the bends 
where the envelope curve can be modeled by a parabola (see the top part of Fig. [l]): 
y n = and z n = with k being the inverse radius of curvature at the bending 

point. We assume that the chain is inextensible that is the distance between neighboring 
sites is constant: \r n +i ~ r n \ = 1 for all n and k. Therefore, when k is not too large 
(|^m — Tn\ > 1 f° r ai l \n — m\ > 1) the excitation dynamics does not depend on the 
curvature of the system in the nearest-neighbor approximation, when J n m = J5 n _ m ,±i. 

In many physical problems, however, the nearest-neighbor approximation is too 
crude. For example, effective dispersive interaction between nonlinear layers in nonlinear 
dielectric superlattices is exponential, exp(— ar), with the inverse radius a depending 
on parameters of the lattice (superlattice spacing, linear refractive index and so on) [21 . 
In the framework of the Peyrard-Bishop model p5, |26| the nonlinear dynamics of the 
DNA molecule can be described by Eq. (Q) with nonlinear term of the opposite sign 
(because of the repulsive anharmonicity) . In this case ip is the complex amplitude of 
the base-pair stretching vibration ||22|| . The base-pairs (AT, GC) in the DNA molecule 
are usually asymmetric (see, e.g., |]J), and therefore the vibration which corresponds 
to the stretching of the hydrogen bonds connecting the two bases, transfers along 
the chain due to the dipole-dipole interaction decaying with the distance r as 1/r 3 . 
Since the dipole moments of base-pairs are almost perpendicular to the molecular axis 
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Figure 2. Dependence of the spectrum of linear excitations on the curvature of the 
chain for a = 2 and J = 6.4. Localized linear bound states emerge above (symmetric 
states) and below (antisymmetric states) the continuum spectrum of extended states. 
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Figure 3. Profiles of the linear bound states from Fig. || for n = 15. 



0, they are almost parallel to each other and the interaction between dipoles is also 
repulsive {J n>m < 0). Thus, after substitutions t — > — t and J„ im — ► — J n<m , this model 
of the DNA molecule will be described by Eq. (0). We have investigated both cases, 
Jn,m = J exp (— a\r n — r m \) and J n>m = J\r n — r m \~ s , with qualitatively the same results. 
In this paper, however, we discuss only the first case in detail. It is known [2^J that in 
the case of straight rod-like chain (« = 0) such NLS model exhibits bistability in the 
spectrum of nonlinear stationary states for a < 1.7. To distinguish these effects from 
the finite curvature effects we use a = 2 in what follows. 

For stationary states, ip n (t) = 4> n (A) e lAt , where A is the nonlinear frequency and 




Figure 4. Number of excitations N(A) for the straight chain (k = ; dot-dashed 
line) and curved chain (re = 3 ; so/id Zmes for stable states, dashed lines for the states 
unstable with respect to the breathing mode, and dotted line for the states unstable 
with respect to the Peierls mode) for a = 2 and J = 6.4. In the inset shapes of the 
symmetric (s) and asymmetric (as) stationary states for N — 3.6 are shown. 



</> n (A) is the excitation amplitude of the n-th site, Eq. (Q) reduces to the system of 
nonlinear algebraic equations 

E J nm(t>m + 0^ = (2 + A) 0„ , (3) 

which can be solved, for instance, by the Newton-Raphson iterations. 

In Fig. ^| we demonstrate that in the linear limit, i.e. when the term 0^ in Eq. (|3|) 
is neglected, the bending creates at least one linear localized state. The profile of the 
first bound state, which is always present, is shown in Fig. |||a. When the curvature 
of the chain increases, new localized states are emerged from the continuum spectrum; 
their profiles are shown in Fig. |3|. As it has been recently proven for another but similar 
system , the number of these bound states can exceed any positive integer provided 
the chain is curved "enough" . 

However, in the present paper we shall study the case of smooth bendings for which 
only the first bound state is present. In Fig. f| we plot the dependence of N(A) obtained 
numerically for both, the straight rod-like and curved chains. As may be inferred from 
this figure, several new features arise as a consequence of the finite curvature. 

Specifically, there appears a gap at small A in which no localized solution can exist. 
As can be seen from Fig. |2], this gap originates from the existence of the linear two-hump 
bound state discussed above. 

In the case of curved chain there are two branches of stationary states: a branch 
of symmetric localized excitations (s), which exists for all values of N, and a branch of 
asymmetric localized excitations (as), which exists only for N > N t h{n)- The threshold 
value N t h(K) decreases when k decreases and vanishes in the limit k — > 0. The symmetric 
stationary state has a two-humped shape (evolved from the linear two-hump localized 
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Figure 5. Energy of symmetric (s) and asymmetric (as) nonlinear excitations versus 
the curvature k for a = 2, J = 6.4, and N = 4. 



mode) while in the asymmetric state the maximum is shifted either to the left or to the 
right from the center of symmetry n = (see inset in Fig. |J). 

In contrast to the case of straight rod-like chain (see dot-dashed line in Fig. |j), the 
dependence N(A) for the symmetric stationary states in curved chain is non-monotonic. 
That is, there exist an interval of N in which three different symmetric states coexist 
for each excitation number; two of them are usually stable [J27], |25fl . Such bistability 
phenomenon exists in the straight rod-like chain only for smaller values of a < 1.7. Thus, 
one may conclude that the bending of the chain facilitates the bistability of nonlinear 
excitations. 

The energy of the symmetric as well as asymmetric localized states is monotonically 
decreasing function of the curvature (see Fig. Thus, one may expect that the 

nonlinear localized excitation may facilitate bending of flexible molecular chain. For 
iV > Nth(ii), when symmetric and asymmetric stationary states coexist, the asymmetric 
state is always energetically more favorable. Taking into account that in the linear limit 
the ground state of the model is symmetric with respect to the center of symmetry n = 
one can conclude that combined action of the finite curvature and nonlinearity provides 
the symmetry breaking in the system. 

Apart from the smooth bending we studied the case of the wedge: y n = \x n \ tan(#) , 
and obtained qualitatively the same results. We checked also that the exact form 
of the dispersive interaction is not of serious concern here. To this end we studied 
the case when the dispersive interaction has the power dependence on the distance: 
J n ,m = J\r n — r m \~ s . Qualitatively the same results were obtained for this case, too. 
The symmetry breaking, which we observe, means that there is a bistability in the 
spectrum of nonlinear excitations: the excitation can be localized either to the left or 
to the right of the bending. By applying an external periodic field one can achieve 
switching between these two states. In this way we may expect a new resonance to 
appear in the response of the system to a periodic driving. 
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3. Analytical approach and minimal model 

To gain a better insight into the symmetry breaking mechanism, let us consider the 
continuum limit of the discrete NLS model given by Eqs. (fj])-(0). We are interested 
here in the case when the characteristic size of the excitation is much larger than the 
lattice spacing. It permits us to replace if) n (t) by the function if)(£,t) of the arclength 
I which is the continuum analogue of n. Using the Euler-Mclaurin summation formula 
we get 



H 



2m(£) 



dip 2 i rr//A| ( |2 1 



di 



+ t/(£)H 2 --H 4 Ut. (4) 



Here m 1 (£) = f (£' — £) 2 Je^>d£' is the inverse effective mass of the excitation, and 

— CO 

CO 

U(£) = — f Je,e'd£' is the energy shift due to the dispersive interaction. Near the 

— CO 

bending point the number of neighbors available for excitation is larger than far away 
where the curvature is small. Near the bending the effective potential, U(£), exhibits a 
double-well shape. This fact together with the resulting double-well spatial dependence 
of the effective mass, m(£), are the reasons why the excitation is localized. The 
attractive potential U(£) consists of two symmetric wells whose positions ( ± a) and 
depths (A = U(\£\ — > oo) — U(a)) are determined by the range of the dispersive 
interaction and the peculiarities of the bending geometry (see Fig. [I]). E.g., in the 
case of the wedge, y n = \x n \ tan(0) , we get for 9 < 7r/4 

U(£) = - J -(2- e-W) - - { e ~ aWC ° Sm WhGn 9 < 1/2 ' 
W a y ' a \ a \£\ #i(a|*|) when 9 -> vr/4 , 

where K n (z) denotes the modified Bessel function The distance between the 

potential minima 2a ~ (1 + 9 2 )/a and their depth A w 4eJ9 2 /a increase when the 
wedge angle 9 —>■ n/4. Note that the Hamiltonian similar to Eq. (f|) with U(£) = was 
introduced in Ref. |3l| to study the effects of the DNA bending on breather trapping. 



The step function dependence of the effective mass m(£) was there assumed. 

The NLS model @ is still very difficult to solve analytically. But it is interesting 
anyway to find a model — chosen as simple as possible — which gives qualitatively the 
same results as the original model. Such a model can be obtained by using step-wise 
approximation 

2 7 

U(£) = A y (H(£ -aa + w/2) - H(£ - aa - w/2)) , (5) 

where H(x) is the Heaviside step-function, w is the width of the potential well. The 
similar expression can be written for the effective mass m{£). In the step- wise approach 
the stationary equation which corresponds to the Hamiltonian (^) can be solved and 
the corresponding eigenvalue problem can be explored. However, results obtained in this 
manner look very cumbersome. Therefore we shall make a further simplification by 
applying the limit w — ► 0, keeping the product wA = e constant in Eq. @. In this 




Figure 6. The dependence N(A) for symmetric (s) and asymmetric (as) nonlinear 
excitations in the two-impurity minimal model (Bh with e — 0.7 and a = 1. 



way we arrive at the following minimal model 

. dip d 2, ip 



+ 



+ e(5(£ -a) + 5{£ + a))ip + |V| 2 V = 



(6) 



dt d£ 2 

where the gauge-transformation ip — > ipexp (—i2Jt/a) was used. Introducing into Eq. 
(^) the stationary state ansatz ip(£,t) = <$(£) exp(iA) where cj){£) — > when \£\ — > oo , 
we find that 



' v / 2Asech(v / A(£ 



'2A(1 -m) 



nd 



A 



V 2 — m 
v / 2Asech( v / A(^-^i) 



in 



when £ < —a, 
(£ - £ ) | m ) when \£\ < a, 
when £ > a 



where pq(w | m) (p, q = c, s, d, n) are the Jacobi elliptic functions with the modulus 
m f50fl. The parameters £±i,£q, and m are determined from the relations 



sech V A (a =f £±i 



1 — m 



nd 



A 



rn 



m 



-(a=f £ \ m) 



F(£ ,m) = \Jl — m (1 — m) sd 2 (w + | m) + msd(w + | m) cn(w + | m) 
\Jl — m(l — m) sd 2 (u^ \ m) + m sd(u_ | m) cn(u_ | m) = 0, 



1 — m (1 — m) sd 2 (w-|- \ m) + m sd(u.|- | m) cn(u_|- | m) 



m 



A 



(7) 



where u± = J A/ (2 — m) (a ± £$). Using the above expressions one can evaluate the 
number of excitations which corresponds to the stationary solution 



N =4 (v 7 ! - e) + 2 



A 



E(am(u a ) | m) 



(8) 



CT=±1 
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Figure 7. Switching from the symmetric state to the asymmetric state at a = 2, 
J = 6.4, k = 4, and N = 3.96. 

where am(w) is the amplitude and E(<f) | m) is the elliptic integral of the second kind 

The equation F(£q, m) = always has a solution £q = 0. This corresponds 
to the symmetric two- hump wave function <f>(£). In addition, £ 7^ solutions 
appear for A > A c with the threshold value A c being determined from the condition 
dF/d£o = for £ = 0. In other words, a symmetry-breaking bifurcation to a doubly 
degenerate branch of asymmetric solutions occurs in the point A = A c . It follows 
from Eq. (0) that the symmetry-breaking bifurcation is determined by the equation 
(1 — 2sn 2 +m sn 4 ) y/1 — 2m sn 2 + m 2 sn 2 = (1 — m) sn cn with sn = sn (^J K/(2 — m) \ m^J . 

The results of the analytical consideration of the nonlinear eigenvalue problem 
based on Eqs. (0) and @ are shown in Fig. [f| As in the original model the stationary 
state (s) is unique and symmetric for small number of quanta N. But when the number 
of quanta exceeds some critical value there are two stationary states: symmetric (s) 
and asymmetric (as), with the asymmetric state being energetically favorable. Thus, 
similarly to the original model, the minimal model demonstrates a symmetry breaking 
effect (localization in one of the wells). It is worth noting that a closely related 
phenomenon was observed in Ref. [|3^] where nonlinear electromagnetic waves in a 
symmetric planar waveguide were studied. But in contrast to our case the nonlinear 
Schrodinger model with a step-wise nonlinearity was considered. 

4. Dynamical properties of excitations 

The dynamical properties of nonlinear excitations include both, the behavior of immobile 
stationary excitations when they lose stability, and a scattering of the moving excitations 
by the bending. 
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Figure 8. Velocity v sc of the moving part of the excitation after the passage of the 
center of the parabolic chain versus the curvature k. Here a = 2, J = 6.4, N = 2.39 
and initial velocity v = 0.37. 



The first issue has been considered in details in Ref. [|33| where we have shown 
that the symmetric two-humped stationary state can lose stability in two ways, due 
to softening of the antisymmetric Peierls internal mode and due to softening of the 
symmetric breathing mode. A typical evolution of the symmetric stationary state in 
the case when it becomes unstable with respect to the Peierls mode is shown in Fig. [5]. 
One can see that the two-humped initial symmetric state evolves into an asymmetric 
state with an excited breathing internal mode. 

In this paper we consider in more details another item which is very intriguing 
from the viewpoint of possible applications: the scattering of moving excitations in 
their passage through the center of the bending. To accomplish it, we first obtained a 
stationary moving excitation by pushing (i.e., taking ip n {f)) = n (A) e lkn with some k 
and A) of the initially unmoving stationary state <fi n found numerically for the straight 
chain as described in Sec. 2. Eventually, evolving over the course of 5000 time units, 
the excitation reached almost stationary form. Substituting it into the curved chain 
far away from the bending, we found out that the destiny of the excitation in passing 
though the bending region is strongly dependent on the curvature k of the chain. 

In Fig. [| we plot the eventual velocity v sc (k) of the excitation with N = 2.39 and 
initial velocity v = 0.37. One can see that for small k (k, < 0.2) the bending does not 
practically change the excitation velocity. However, for 0.2 < k < 1.07 the transmission 
becomes clearly inelastic, leading to decreasing of the excitation velocity. The excitation 
passes the bending region leaving behind a part of its energy for exciting of the linear 
bound state (see Fig. |^a). 

In the interval 1.07 < k < 2.1 the excitation is captured by the bending (see 
Fig. |9|b), so that in this case the bending clearly manifests itself as an effective trap. 



Similar phenomenon has been previously considered [31] as an effective breather trapping 
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Figure 9. Density plot of |^ n (f)| 2 which demonstrates scattering of the moving 
excitation from Fig. || for different values of k: (a) 1.05, (b) 1.5, (c) 2.5, and (d) 
3.0. 



mechanism for DNA transcription in the case of curved DNA. However, in Ref. [31] the 
bending has been accounted effectively, postulating some modification of the coupling 
constants around it. 

Finally, for 2.1 < k < 2.6 the excitation is reflected from the bending, however a 
significant part of the initial excitation energy remains captured near the bending region 
(see Fig. |9]c). In contrast, for k > 2.6 the main part of the excitation energy is reflected 
from the bending and only small part is captured (see Fig. |9]d). 



5. Conclusions 



In conclusion, we have shown that the interplay of curvature and nonlinearity in the 
chain of nonlinear oscillators can lead to the qualitative effects, such as symmetry 
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breaking of the nonlinear excitations and their trapping by the bending. We have 
demonstrated that the energy of nonlinear excitations decreases with the increasing of 
curvature of the chain and thus one can expect that the presence of nonlinear localized 
excitation in the chain may facilitate its bending, if only the chain is flexible enough. 

These effects are tolerant to the choice of the coupling coefficients J nm . In 
particular, qualitatively the same properties were obtained for the case of the dipole- 



dipole dispersive interaction with J, 
Davydov-Scott model | 



r„ - r r . 



Specifically, in the case of the 



19| with the dipole-dipole coupling between sites |20f for the 



parameters which characterize the motion of the amide-I excitation in proteins (they 
correspond to N = 0.64 in our model) decreasing of the ground state energy due to a 
trapping by the bending with k — 1 (k — 3) comprise 10% (46%) of the whole excitation 
energy. This may be also important for the dynamics of DNA molecule and in particular 
for DNA melting where in accordance with Peyrard-Bishop model [25| |26| the nonlinear 
properties of base-pairs vibrations play a crucial role. 



Acknowledgments 

We thank A.C. Scott and M. Peyrard for the helpful discussions. Yu.B.G. and S.F.M. 
would like to express their thanks to the Department of Mathematical Modelling, 
Technical University of Denmark where the major part of this work was done. S.F.M. 
acknowledges support from the European Commission RTN project LOCNET (HPRN- 
CT- 1999-00163). 



References 

Saenger W 1984 Principles of Nuclear Acid Structure (Berlin: Springer- Verlag) 
Iijima S 1991 Nature 354 56 

Bar-Ziv R and Moses E 1994 Phys. Rev. Lett. 73 1392 
Nelson P, Powers T and Seifert U 1995 Phys. Rev. Lett. 74 3384 

Meunier J, Langevin D and Boccara N (Eds.) 1987 Physics of Amphiphilic Layers (Berlin: Springer- 
Verlag) 

Imry Y 1997 Introduction to Mesoscopic Physics (New York: Oxford University Press) 
Soukoulis C M (Ed.) 1993 Photonic Band Gaps and Localization (New York: Plenum Press) 
Reiss C 1995 Nonlinear Excitations in Biomolecules Ed.: Peyrard M (Berlin: Springer, Les Ullis: 

Les Editions dc Physique) p 29 
Hcumann H, Ricchetti M and Werel W 1988 EMBO J. 7 4379 
Burley S K, 1994 Nature Struct. Biology 1 638 
Kim J L, Nikolov D B and Burley S K 1993 Nature 365 520 
Hamada N, Sawada S and Oshiyama A 1992 Phys. Rev. Lett. 68 1579 
Menon M and Srivastava D 1997 Phys. Rev. Lett. 79 4453 
Lenz G, Talanina I and de Sterke C M 1999 Phys. Rev. Lett. 83 963 
Kirczenow G 1989 Phys. Rev. B 39 10452 

Gaididei Yu B, Malysheva L I and Onipko A I 1992 J. Phys.: Condens. Matter 4 7103 
Dandoloff R, Villain- Guillot S, Saxena A and Bishop A R 1995 Phys. Rev. Lett. 74 813 
Saxena A and Dandoloff R 1997 Phys. Rev. B 55 11049 
[16] Saxena A, Dandoloff R and Lookman T 1998 Physica A 261, 13 



Effects of finite curvature on soliton dynamics 



13 



Harayama T, Davis P, Ikeda K S 1999 Phys. Rev. Lett. 82, 3803 

Davydov A S 1985 Solitons in Molecular Systems (Dordrecht: Reidel Publishing Company) 
Scott A C 1992 Phys. Rep. 217 1 
Cruzeiro-Hansson L 1998 Phys. Lett. A 249 465 

Yakushevich L V 1998 Nonlinear physics of DNA (New York: John Wiley & Sons) 

Gaididei Yu B, Mingaleev S F, Christiansen P L, Johansson M and Rasmussen K 1998 Nonlinear 

Cooperative Phenomena in Biological Systems Ed.: Mattson L (Singapore: World Scientific) 

p 176 

Golestanian R, Liverpool T B 2000 Phys. Rev. E 62 5488 

Gaididei Yu B, Christiansen P L, Rasmussen K and Johansson M 1997 Phys. Rev. B 55 R13365 

M. Peyrard and A. R. Bishop, Phys. Rev. Lett. 62, 2755 (1989). 

T. Dauxois, M. Peyrard and A. R. Bishop, Phys. Rev. E 47, 684 (1993). 

Johansson M, Gaididei Yu B, Christiansen P L and Rasmussen K 1998 Phys. Rev. E 57 4739 
Exner P 2000 Bound states of infinite curved polymer chains arXiv:math-ph/0010046] 



Gaididei Yu B, Mingaleev S F, Christiansen P L and Rasmussen K 1997 Phys. Rev. E 55 6141 
Abramowitz M and Stegun I 1972 Handbook of Mathematical Functions (New York: Dover 
Publications) 

Ting J J-L and Peyrard M 1996 Phys. Rev. E 53 1011 

Akhmediev N N 1982 Sov. Phys. JETP 56 299 (1982 Zh. Eksp. Teor. Fiz. 83 545). 
Gaididei Yu B, Mingaleev S F and Christiansen P L 2000 Phys. Rev. E 62 R53 
Mingaleev S F, Christiansen P L, Gaididei Yu B, Johansson M and Rasmussen K 1999 J. Biol. 
Phys. 25 41 



